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Summary 


The  first  collection  of  reports  written  by  L.C.  Woods  for  the  UT  Institute  for 
Advanced  Technology,  University  of  Texas  at  Austin,  contained: 


1.  The  Current  Melt- Wave  Model,  October,  1995. 

2.  The  Contact-Spot  Model  of  Transition,  December,  1995. 

3.  Paper  for  the  8th  EML  Symposium,  January,  1996. 

4.  An  Armature/Rail  Instability,  February,  1996. 

5.  Three-Dimensional  Effects  on  the  Melt- Wave  Model,  May,  1996. 

6.  Boundary  Conditions  Across  a  Rail/Armature  Contact  Surface, 
June,  1996. 


The  present  collection  adds  four  more  reports  to  the  list.  They  are  numbered 
as  in  the  actual  papers,  the  number  7  being  assigned  to  the  first  collection.  The 
addition  reports  are: 

8.  Viscous  Effects  on  the  Liquid  Interface  Between  Armature  and  Rail, 

November,  1996. 

9.  The  Distribution  of  the  Electromagnetic  Force,  November,  1996. 

10.  Testing  the  Melt- Wave  Model,  January,  1997. 

11.  Investigation  of  Transverse  Striations  on  Aluminium  Rails, 
(Co-author,  Dr.  Andrew  Yeoh),  November,  1996.  This  paper  is  to  be  pre¬ 
sented  at  the  6th  European  Symposium  on  EML  Technology,  to  be  held  at 
The  Hague  in  May,  1997. 

A  brief  survey  of  their  contents  and  relationships  follow. 

8.  This  report  had  several  objectives,  the  first  of  which  was  to  determine 
the  rate  at  which  the  armature  is  melted  by  viscosity  and  to  use  this  to  find  the 
thickness  of  the  layer  of  aluminium  deposited  on  the  copper  rail  in  circumstances 
when  ohmic  heating  is  not  important.  To  this  end  it  was  necessary  to  determine 
the  velocity  distribution  in  the  liquid  layer  at  the  armature/rail  interface  and 
also  to  estimate  the  magnitude  of  the  pressure  gradient  in  the  layer  due  to  the 
applied  electromagnetic  force.  This  latter  problem  was  not  straightforward  and 
resulted  in  the  work  to  be  reported  in  9.  A  second  objective  was  to  investigate  the 
stability  of  the  liquid  layer  under  the  viscous  and  electromagnetic  forces.  Were 
it  found  to  be  marginally  unstable,  it  was  conjectured  that  this  might  provide 
an  explanation  for  oscillations  that  have  been  observed  in  the  muzzle  voltage  in 
some  lAT  exj)eriments.  But  the  frequency  of  the  oscillations  predicted  by  theory 
proved  to  be  much  too  small  to  account  for  the  observations. 

A  third  aim  was  to  derive  an  expression  for  the  solidification  time  4,  i.e.  the 
time  that  il,  taluks  the  melted  layei'  to  solidify  on  the  rail,  4  b6ing  measured  from 
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the  instant  that  the  melt  leaves  the  rear  of  the  armature.  For  thin  layers,  vis¬ 
cous  heating  is  important  since  it  tends  to  reheat  the  layer.  This  work  was  in 
anticipation  of  the  investigation  of  the  instability  to  be  described  below  in  11. 
Finally  ohmic  heating  was  restored  and  a  comparison  made  between  between  the 
two  heating  mechanisms.  It  was  found  that  for  liquid-layer  thicknesses  exceed¬ 
ing  10  /rm,  ohmic  heating  is  typically  more  important  than  viscous  heating,  the 
difference  resulting  from  the  fact  that  ohmic  heating  is  very  localized. 

9.  To  determine  the  force  acting  on  the  liquefied  armature  material  at  the 
rail/armature  interface,  it  is  necessary  to  know  how  the  total  electromagnetic 
force,  F  =  \L'P,  where  I  is  the  total  electric  current  and  V  is  the  gradient  of 
the  self-inductance,  is  distributed  over  the  rear  of  the  armature.  In  particular  the 
magnitude  of  the  force  at  the  rear  corners  of  the  interface  liquid  is  required.  If 
the  force  near  these  corners  is  rather  larger  than  the  average,  it  would  tend  to 
accelerate  the  liquid  layer  forward,  ahead  of  the  solid  part  of  the  armature.  The 
Biot-Savart-Lorentz  (BSL)  formula  was  used  to  deduce  an  analytic  expression  for 
this  force  and  revealed  that  there  is  only  a  small  variation  in  its  value  over  the 
rear  of  the  armature,  30%  being  typical. 

A  convenient  analytic  expression  for  U  as  a  function  oi  n  =  hr/wr  and  s  = 
t/wr,  where  F  is  the  distance  between  rails  of  width  Wr  and  thickness  t.  This  is 
not  exact,  since  it  was  assumed  that  the  current  distribution  was  uniform.  But 
for  small  values  of  5,  it  appears  to  be  a  good  approximation,  and  in  agreement 
with  observation.  Finally  it  was  shown  that  the  claim  that  the  original  Amperian 
law  for  the  force  between  two  current  elements  (which  is  directed  along  the  vector 
joining  the  elements),  gives  ‘very  nearly  the  same  acceleration  force’  as  the  BSL 
law  was  investigated.  Typically,  the  Amperian  law  overestimates  this  force  by 
50%,  .so  that,  apart  from  lacking  a  physical  mechanism,  this  law  has  no  merit 
from  an  experimental  point  of  view. 

10.  This  report  returns  to  the  melt-wave  model  described  in  earlier  reports 
and  is  concerned  with  the  question  ‘what  observations  could  be  made  to  test 
it?’  The  difficulties  of  measurement,  at  least  in  present  experiments,  restrict  the 
observations  to  the  thickness  6m  of  the  aluminium  layer  deposited  on  the  copper 
rails  and  to  the  depth  h  of  the  aluminium  layer  removed  by  melting  from  the 
armature.  As  these  two  lengths  are  closely  related  by  continuity,  and  depend  on 
the  ratio  of  the  melt  wave  velocity,  u„,,  to  the  armature  velocity,  F,  we  could 
using  them  to  deduce  values  of  w,„,  but  direct  measurement  of  this  speed  would 
be  difficult.  At  this  stage,  we  are  left  with  6m  as  providing  the  only  link  between 
theory  and  experiment. 

In  report  1,  the  gap  depth  h  was  determined  from  the  diffusivity  of  the  electric 
current  on  the  assumption  that  the  mechanism  was  essentially  adiabatic.  In  2 
doubt  s  about  the  validity  of  this  were  advanced,  it  being  argued  that  the  thermal 
diffusivity  should  have  the  dominant  role.  The  variation  of  6m  with  distance  x 
along  the  rail  follows  quite  different  laws,  depending  on  whether  it  is  thermal 
or  ohmic  diffusivity  that  determines  the  value  of  h.  One  lAT  experiment  clearly 
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favoured  thermal  diffusivity,  leaving  the  difficult  task  of  determining  the  constant 
of  proportionality.  In  the  report  a  semi-empirical  treatment  of  melting  led  to  an 
expression  for  dm  in  close  agreement  with  the  experiment.  More  experiments  in 
which  6m  is  measured,  are  required  to  support  to  the  theory. 

11.  This  paper  (written  jointly  with  Dr.  Andrew  Yeoh  of  lAT)  is  to  be 
presented  at  the  6th  European  Symposium  on  EML  Technology,  in  May,  1997. 
It  deals  with  a  curious  phenomenon  that  occurs  when  both  the  rails  and  the 
armature  are  made  of  aluminium.  In  certain  conditions,  periodic,  transverse 
striations  are  found  in  the  rail  surface,  with  a  wave  length  varying  between  0.5 
mm  and  about  5mm.  It  appears  that  the  melted  aluminium  at  the  contact 
surface  has  solidified  into  waves  with  their  crests  at  right  angles  to  the  motion  of 
the  armature.  It  is  likely  that  these  waves  are  generated  by  an  instability  in  the 
sheared  fluid  layer  separating  the  rail  and  armature,  at  a  point  just  to  the  rear 
of  the  armature.  The  report  gives  an  account  of  the  experimental  conditions  and 
of  the  type  of  instability  that  may  be  involved  in  the  phenomenon.  The  work  on 
the  solidification  time,  reported  in  8,  was  adapted  to  show  that  the  time  scale  of 
the  growth  of  the  instability  is  less  than  that  for  solidification,  which  means  that 
the  waves  have  time  to  develop  before  the  liquid  state  is  lost. 
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(Second  Collection  of  Progress  Deports  on  Desearch  into 
Electromagnetic  Launchers 


Viscous  Effects  on  the  Liquid  Interface 
Between  Armature  and  Dail 


November  1996 


Professor  L.  C.  Hoods 
University  of  Oxford 
Mathematical  Institute 


Summary 


In  some  circumstances  viscous  heating  at  the  armature/rail  interface  is  more 
important  than  ohmic  heating.  In  this  report  is  investigated  the  influence  of 
viscosity  on  armature  melting  in  the  absence  of  ohmic  heating.  The  rate  at 
which  the  armature  is  ‘consummed’  by  viscosity  is  determined  and  used  to  find 
a  value  for  the  thickness  of  the  layer  of  aluminium  deposited  on  the  copper  rail. 
The  calculation  requires  us  to  determine  the  velocity  distribution  in  the  liquid 
layer  at  the  interface  and  also  to  estimate  the  magnitude  of  the  pressure  gradient 
in  the  layer  due  to  the  electromagnetic  force. 

It  is  shown  that  there  is  a  length  ^  for  the  liquid  layer,  measured  in  the 
direction  of  motion,  at  which  equilibrium  is  achieved.  If  the  layer  is  shorter  than 
this,  the  EM  forces  accelerate  it  forwards  and  if  it  is  longer,  it  falls  back  to  its 
equilibrium  position.  There  is  a  natural  frequency  uj  for  oscillations  in  £.  It  was 
initially  hoped  that  this  frequency  would  provide  an  explanation  for  oscillations 
that  have  been  observed  in  the  muzzle  voltage  in  the  lAT  experiment,  LeOAT 
33,  but  u  is  much  too  small.  It  may  have  other  applications. 

The  report  also  contains  a  theory  for  the  time  ^  that  it  takes  the  melted 
layer  to  solidify  on  the  rail,  where  the  time  is  measured  from  the  instant  that  the 
layer  leaves  the  rear  of  the  armature.  For  thin  layers,  viscous  heating  becomes 
important  since  it  tends  to  reheat  the  layer.  The  significance  of  this  result  is  that 
for  instabilities  in  the  layer  downstream  of  the  armature  to  be  visible  on  the  cold 
rails  after  an  experiment,  they  must  have  a  growth  rate  smaller  than  4. 

Finally  ohmic  heating  is  restored  and  a  comparison  is  made  between  between 
the  two  heating  mechanisms.  It  is  shown  that,  provided  the  thickness  of  the  liquid 
layer  exceeds  a  critical  value,  about  10 /xm  for  typical  values  in  lAT  experiments, 
ohmic  heating  is  more  important  than  viscous  heating.  It  is  explained  why  this 
result  differs  from  a  conclusion  recently  presented  at  lAT  by  Dr.  D  Hildenbrand. 


1.  Introduction 

An  example  of  purely  viscous  melting  occurs  in  experiments  conducted  at 
lAT  aimed  at  determined  the  gouging  threshold.  The  armatures  used  have  three 
sections,  the  central  one  being  electrically  isolated,  and  therefore  heated  by  vis¬ 
cosity  alone.  Our  objective  here  is  to  determine  the  rate  at  which  viscous  heating 
occurs,  with  allowance  made  for  thermal  losses  to  the  rails.  An  expression  is 
derived  for  the  thickness  of  the  layer  deposited  on  the  rails. 

Figure  1  illustrates  the  two-dimensional  model  that  shall  be  adopted.  The 
armature  is  pressed  on  to  the  moving  rail  by  a  force  F,  and  is  melted  by  the 
heat  generated  by  friction  in  the  liquid  layer.  This  layer  has  a  thickness  4  and 
an  average  velocity  when  it  solidifies  on  the  rail,  its  thickness  is  <5^  and  its 
velocity  is  V.  If  Vm  is  the  velocity  of  the  solid  armature  towards  the  melt  front 
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(1) 


and  a  is  the  width  of  the  armature  feeding  into  this  front,  by  continuity 

In  the  following,  we  shall  use  the  steady-state  energy  equation  to  determine 
Vm  and  to  complete  the  calculation  we  need  an  expression  for  v'.  The  velocity 
distribution  in  the  liquid  layer  is  determined  from  the  equation  of  fluid  motion, 
which  in  turn  requires  us  to  evaluate  the  EM  force  acting  on  the  layer. 

2.  Velocity  of  the  Melt  Front 

We  shall  start  with  the  energy  equation 

PCpi^  +Vm^)T  -h  T]f  +  puv'^  -  pVmLm^{6s),  (2) 

ay' 

where  p  is  the  density  of  the  liquid  aluminium,  u  is  its  kinematic  viscosity,  K  is 
the  thermal  conductivity,  r]  is  the  electrical  resistivity,  j  is  the  current  density, 
Lm  is  the  enthalpy  of  fusion,  v'  =  dv/dy,  and  6{ym)  is  the  delta  function, 
y  =  6s  being  the  location  of  the  melt  front  (see  Fig.  1).  At  first  we  shall  omit 
the  ohmic  heating  term;  we  shall  also  ignore  heating  losses  in  the  OV-direction, 
reducing  (2)  to 

pCpVm^  =  +  pl'v’'^  -  pVjnLm^{8s),  (3) 

in  steady  conditions.  The  next  step  is  to  integrate  this  expression  over  the  in¬ 
terface  surface,  0  <  a;  <  a,  where  x  is  measured  from  the  leading  edge  of  the 
armature  in  electrical  contact  with  the  rail.  Let  denote  the  melting  tempera¬ 
ture  of  the  liquid  aluminium,  and  To  the  temperature  on  the  outside  of  a  rail  of 
thickness  A^,  then  ignoring  ohmic  heating  in  the  rail,  its  temperature  is  given  by 

T-ro  =  (l-y/A,)(T,„-To)  (-A,  <  y  <  0).  (4) 

where  we  have  assumed  that  the  liquid/rail  surface  is  on  y  =  0. 


I  solid 
V  p  armature 


Rail - ►  V 


Figure  1:  Armature  melted  by  a  moving  rail 
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Figure  2;  Velocity  profile 


Since  is  zero  outside  4  <  y  <  oo,  u  is  zero  outside  0  <  y  <  ^  and  Kr  is 
zero  outside  <  y  <0,  the  integral  of  (3)  over  — A^  <  y  <  oo  gives 

PCpVmiTm  -  To)  =  -  Tq)  +  pu  dy  —  pVmLm-  (5) 

Let  Q  =  c(r„  -  To)  +  i„,  Hs  -^(T„  -  To), 

pQK 

then  (5)  can  be  written 

^  ,2 

'’”’  =  q/o  ”  («) 

To  continue,  we  need  an  expression  for  the  velocity  gradient  in  the  liquid  layer. 

3.  The  Velocity  Distribution  in  the  Liquid  Layer 

We  shall  chose  a  reference  frame  for  the  velocities  in  which  the  velocity  of  the 
solid  part  of  the  armature  is  stationary  and  therefore  the  rail  is  moving  along  the 
positive  OV-axis  with  speed  V.  The  equation  of  fluid  motion  in  this  frame  is 

dv  /  \ 

W  (7) 

where  p  is  the  fluid  pressure  and  8^2^,  is  the  magnetic  pressure.  In  the  steady 
state,  with  motion  along  OX.  ^  this  equation  gives 


P-  ^  ^ 

dy^  ’  + 


(8) 


With  the  boundary  conditions  i;  =  V  at  y  =  0  and  n  =  0  at  w  =  &  f8)  has 
the  solution  ’ 

^  ~  2^y(y  ~  4)  +  v(i  —  y/6s). 


Later  we  shall  require  the  gradients 


-  v/6s ,  u;  =  -iP4  -  \//4  ^ 


(9) 

(10) 
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where  the  subscript  ‘a’  denotes  the  solid  armature  surface,  y  =  4?  and  ‘r’  the  rail 
surface,  ?/  =  0.  The  average  velocity  is 


Vs  =  \v-  (11) 


The  minimum  velocity  occurs  at  v'  =  0,  i.e.  at  y  =  yo,  where 

!/.=  |4  +  ^.  (12) 

This  is  above  the  midplane,  as  indicated  in  Fig.  2.  The  velocity  at  this  point  is 


(13) 

Also 

Hence  from  (6) 

(14) 

It  follows  from  (1),  (11),  and  (14)  that  8s  is  determined  by  the  quartic 

^P(l  -  ^P)8t  +  +  aH6s  -  =  0,  {£  =  ua/Q) ,  (15) 

which  has  just  one  positive  root.  When  4  is  evaluated,  6m  follows  from 

=  =  (16) 


4.  The  Electromagnetic  Force 

The  electromagnetic  force  on  the  rear  of  the  armature  varies  slowly  from  a 
minimum  value  at  the  centre  of  the  armature  to  a  maximum  at  the  rail,  but  it  is 
a  slow  change,  with  a  maximum  typically  only  about  30%  or  so  of  the  minimum 
value  [1].  In  terms  of  the  inductance  gradient  L'  and  the  railgun  current  I,  the 
total  force  is  P  =  \L'P.  This  is  due  to  the  difference  in  the  magnetic  pressure 
across  the  armature.  We  shall  distribute  this  uniformly  over  the  armature,  ig¬ 
noring  the  small  variation  Just  described.  Let  6  denote  the  distance  between  rails 
of  width  w.  Then  the  fraction  8s /b  of  F  applies  to  each  of  two  liquid  layers  of 
surface  area  8sW.  The  situation  is  correctly  represented  if  we  assume  that  there 
is  a  pressure  (8sF /b) / {8sw)  at  the  rear  of  the  liquid  layer  and  zero  pressure  at 
the  front,  since  we  are  concerned  only  with  the  pressure  difference.  The  gradient 
follows  on  dividing  by  a  (see  Fig.  1).  Hence  the  magnetic  pressure  gradient  in  the 
conducting  liquid  is  F/{bwa)  directed  along  —OX.  We  shall  ignore  the  material 


10 


pressure  p  since  this  is  unlikely  to  deviate  much  from  atmospheric.  Hence  the 
parameter  P  defined  in  (8)  is 


P  = 


F 

2pubwa 


(17) 


Also  from  (8) 


d  1 

rBK 

F 

dx' 

bwa 

Let  us  assume  that  B  varies  from  to  —Bm  over  the  length  a  of  the  liquid 
layer,  then  an  estimate  of  the  maximum  strength  of  the  magnetic  field  follows 
from 


Bm  = 


///oFU  _ 

2bw  ' 

(18) 


To  make  numerical  estimates  here  and  below,  we  shall  adopt  the  values  (in  SI 
units,  except  where  indicated) 


L'  -  0.47x10  ®,  /  =  4x10®,  a  =  w  =  2  cm,  6  =  2.4  cm,  Ar  =  6.35  mm, 

Lm  =  3.97x10®,  Q  =  9.67x10®,  Tm-To  =  635,  p  =  2410,  i/  =  5.7x10-^ 

Kr  =  402 . 

from  which  it  follows  that  F  =  3.76x10^,  P  =  1.37x10^^  =  0  1096  and 

B,n  =  10T. 

The  value  that  we  expect  for  4  is  ~  10“®,  or  less.  Hence,  the  first  two  terms 
in  (15)  are  much  smaller  than  the  others,  so  that  the  equation  reduces  to 


^2  2aH 

(19) 

Also  2aH/V  <Cl,  so  that 

4  = 

(  ^  «  1.18xlO"*‘(aV')i  . 

(20) 

To  suificient  accuracy,  (16) 

gives  6m  =  |4)  and  therefore 

=  0.59(oK)^  pm. 

(21) 

For  example,  with  a  =  2cm,  1/  =  800ms“‘,  we  get  6m  =  2Apm. 
5.  The  Solidification  Time 


Downstream  of  the  armature,  the  fluid  shear  is  eliminated  by  viscosity  on  a 
time  scale  an  effect  that  can  be  accommodated  by  replacing  (9)  by 

u  =  F(l-|-e-‘/‘^), 
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where  t  is  the  time  that  has  elapsed  since  the  passage  of  the  armature.  Hence  at 
a  distance  x  —  a  =  Vt  from  the  rear  of  the  armature,  the  velocity  is 

(22) 

With  the  typical  value,  v  =  5.7x10“^  m^s“^  for  molten  aluminium,  a  layer 
thickness  of  4  =  5 /im  and  a  rail  velocity  of  800ms“\  we  find  4  ~  44 /iS, 
corresponding  to  a  displacement  oi  Vt„  =  3.5  cm.  We  shall  show  below  that 
the  time  for  the  solidification  of  the  layer  is  longer,  so  in  the  region  between 
the  rear  of  the  armature  and  the  beginning  of  solidification,  viscosity  cannot  be 
neglected.  But  with  a  thicker  layer,  say  4  =  20 /im,  Vti,  is  56  cm  and  many  times 
longer  than  the  ‘solidification’  distance  D  defined  in  Fig.  3.  In  this  case  viscous 
heating  is  negligible. 

In  Fig.  3  we  show  an  armature  extending  over  a  contact  distance  a,  measured 
in  the  direction  of  motion.  The  mixture  remains  liquid  on  the  rail  surface,  until 
a  distance  (a  +  D)  measured  from  the  leading  edge  of  the  armature.  The  rail 
temperature  is  given  by  (4).  In  the  range  a  <  x  <  oo,  the  aluminium  layer  is 
not  ohmically  heated.  In  (2)  the  enthalpy  of  fusion  represents  a  loss  of  energy 
during  melting,  but  in  the  present  application  its  role  is  reversed;  it  is  a  store  of 
energy  to  be  removed  by  solidification.  Hence,  in  steady  conditions,  downstream 
of  the  armature,  (2)  becomes 

Q  ==  Kr-^  +  puv'^  +  fwL^6{x*)  (x*  =  a  +  D),  (23) 

where  from  (22),  v  is  the  average 

u  =  V(1  -  (24) 

Substituting  this  into  (23)  and  integrating  across  the  liquid  layer,  we  get 

At  y  =  0,  dT/dy  =  — (T^  —  To)/Ar;  at  y  =  4)  dT/dy  «  0,  since  we  shall  ignore 
losses  to  the  air.  Hence 

-  To)  =  +  pv{l  - 

Next  we  integrate  from  x  =  a  to  x  =  i*.  It  is  convenient  to  introduce  the 
non-dimensional  distance  2  =  D/Vt^,.  We  find 


Figure  3:  Solidification  Distance,  D 


where 


_  __pvLrn__/^\  q7^-,q-3^ 

^  “  O.kJT^  -  Tni  V  L  )  /)„  ’ 


fXr  > 

2Kr{Tm-TQ)^6s 


and  we  have  used  the  values  given  in  §4. 

With  (25)  solved  for  z,  the  solidification  time  follows  from 

4  =  =  1.88x10= AA{(i  -  je-")  +  ^(1  -  e-=')} 


(26) 


Consider  the  case  4  =  5/im,  Xr  —  6.35  mm,  and  V  =  800  ms“^;  (25)  gives 
z  w  3.55,  since  the  exponential  terms  are  quite  small.  Then  4  =  156  and 
D  =  12.5  cm. 


6.  Acceleration  of  the  Liquid  Layer 

Let  M  —  2m  denote  the  mass  of  the  solid  armature,  m  the  mass  of  each  layer 
and  A  —  saw  the  surface  area  over  which  the  viscous  forces  act  (see  Fig.  4).  The 
viscous  force  on  the  liquid  layer  is  /9^'(t;^  —  v'j.)A  and  that  on  the  solid  armature 
is  —2puv'g^A.  It  follows  that  the  equation  of  motion  of  each  layer — considered  as 
a  whole — is 

=  +  (27) 

and  for  the  solid  armature 

(M-2m){^l  =  -^=^F-2p,v',A.  (28) 

the  subscripts  being  7'  for  liquid  and  ‘s'  for  solid.  Notice  that  twice  (27)  plus 
(28)  gives 

M{^)=-F-2puv',A,  (29) 

provided  that  the  accelerations  are  the  same. 

By  (10),  pi'iv'g^  -  v!r)A  =  puPbsA  =  SsF/2b,  whence  (27)gives 
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Figure  4:  Idealized  Liquid  Interface 


Similarly  we  find 


{M  -  2m)  =  -F  +  {Z6s/b)F  +  2pvsawv/6s . 


'dt 


Hence 


(dv^ 

rdVs 

^dt's 

(31) 

(32) 


where 


G  = 


F  ^  ZSgF  ^  2vmv  ^  8sF 
M  —  2m  b{M  —  2m)  b6s{M  —  2m)  2bm 


Since  m/M  1  and  4/^  L  this  is  approximately 


F  fM6s  ,  2vmv 
~  bSsM  ■ 


(33) 


For  the  solid  and  liquid  parts  of  the  armature  to  have  the  same  acceleration, 
we  require  G  =  0.  Since  both  the  terms  on  the  left-hand  side  of  (32)  are  negative, 
G  >  0  requires  the  liquid  to  accelerate  faster  than  the  solid.  Except  exactly  at 
2bm  =  M6s,  the  second  term  in  (33)  is  much  smaller  than  the  first.  To  verify 
this,  we  compare  F  =  \L'P  with  2umv/b6s,  using  the  typical  values  given  in  §4. 
Hence  it  is  sufficiently  accurate  to  write 


M^2bm  '  M^2pb'wsa 


even  in  the  neighbourhood  of  G  =  0.  It  follows  that 


if 


2pbwsa 


<  ^  liquid  overtakes  solid 
>  '  solid  overtakes  liquid 


(34) 


(35) 
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This  result  is  independent  of  the  armature  geometry.  In  the  case  of  a  parallelep¬ 
iped,  M  =  pbwa,  and  (35)  becomes 


<  1  liquid  overtakes  solid 
>  2  • 

solid  overtakes  liquid 


(36) 


In  (34)  let  m  —  mo,  where  M6s  —  26mo,  i.e.  mo  is  the  equilibrium  mass  of 
the  layer.  Write  m  -  mo  +  x{p6sw),  then  x  is  the  displacement  of  the  liquid  front 
beyond  equilibrium.  Substituting  for  m  in  (34)  and  assuming  x  to  be  small,  we 
get 

„  2Fpwb 
- 

M2 

We  need  to  relate  the  acceleration  of  the  whole  mass  of  the  layer  to  that  of  its 
end,  viz.  x.  For  the  velocity,  we  shall  assume  that  the  additional  mass  required 
is  supplied  uniformly  along  the  whole  layer,  so  that  the  velocity  of  the  centre  of 
mass  is  The  additional  mass  affects  the  acceleration  in  a  similar  manner, 
allowing  the  left  hand  side  of  (32)  to  be  written  x/4.  With  this  assumption,  the 
equation  can  now  be  written 

2  k 

X=‘-Uj\  {u  =  —{2pwbFf).  (37) 


7*  Comparison  with  Experiment 

experiment  [2],  shows  oscillations  in  the  muzzle  voltage  that  could  be 
due  to  modifications  in  the  electric  circuit  due  to  periodic  movement  of  the  liquid 
layer.  There  are  fast  oscillations  in  the  voltage  trace  at  a  frequency  of  /  =  46  kH, 
corresponding  to  an  angular  frequency  of  cj  =  289  kH.  These  occur  where  /  = 
162x103  A,  so  that,  with  L'  =  0.47x10-6,  F  =  6.16x103.  The  mass  of  the  armature 
IS  24.37  grm,  then  with  w  and  b  equal  to  2  cm,  (37)  gives  cj  =  8.9  kH,  a  factor  of  32 
in  error.  We  must  conclude  that  the  hypothesis  concerning  periodic  movements 
of  the  liquid  layer  being  responsible  for  the  observed  voltage  fluctuations  is  not 
correct.  The  inertia  of  the  layer  is  too  large  to  fit  the  observations  and  we  must 
look  elsewhere  for  an  explanation  of  the  phenomenon. 

8.  Ohmic  versus  Viscous  Heating 

Finally  we  shall  determine  an  expression  for  the  ratio  of  the  ohmic  to  the 
viscous  heating.  The  terms  to  be  compared  are  pf  and  pvv'^,  but  in  making  this 
comparison,  we  must  remember  the  different  ranges  over  which  these  processes 
operate.  Viscosity  acts  over  the  whole  of  the  layer,  which  in  the  notation  of  Fig.  3, 
is  the  distance  a.  On  the  other  hand  ohmic  dissipation  tends  to  be  localized  near 
the  trailing  edge  of  the  armature,  even  if  the  conductivity  is  uniform  along  the 
surface.  In  [3]  we  showed  that  the  current  layer  wiis  confined  to  the  very  small 
thickness,  8r  —  77/ p^\  ,  so  that  if  w  is  the  width  of  the  armature,  the  current 
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density  is  related  to  the  total  current  /  by  ji  =  I/SrW.  We  shall  distribute  the 
energy  rjj'^  over  the  whole  of  the  armature  surface  by  multiplying  it  by  6r/a,  so 
the  comparison  is  now  between  {Sr/a)r}j^  and  puv'^.  For  v'  we  write  V/Ss-  Hence 

-jl  =  Qhmic  ^  (3g) 

viscous  \6raw‘^  J  ypvV"^  J  pu  \  w  J  aV 
Using  the  values  cited  in  §4,  we  get 

71  =  9.15x10-^  f—)  4?- 
\wj  aV 

Consider,  for  example,  the  case  w  =  a  =  2 cm,  with  V  =  800 ms  I  =  4x10^  A. 
We  get  TZ  =  2.29x10^“  ^  so  that  7?.  >  1  if  (^  >  6.6  ^m.  At  a  typical  thickness  of 
15  /xm,  TZ  =  5.2,  i.e.  ohmic  heating  is  five  times  greater  than  viscous  heating. 

Turning  to  Dr.  Hildenbrand’s  report  [4],  he  takes  the  energy  ‘sourced’  at  an 
interface  of  area  A  to  be  RP  +  pvv'vA,  where  the  resistance  is  i?  =  r]6slA. 
Therefore  his  version  of  the  ratio  defined  in  (38)  is 

^  pv  \amV )  a  ' 

i.e.  for  the  particular  example  mentioned  above,  he  gives  viscosity  188  times 
greater  importance  than  is  actually  the  case.  His  error  is  to  forget  that,  unlike 
viscous  dissipation,  ohmic  dissipation  is  localized. 
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Summary 


In  order  to  determine  the  force  acting  on  the  liquefied  armature  material  at  the 
rail/ armature  interface,  we  need  to  know  the  distribution  of  the  electromagnetic 
force  accelerating  the  armature  along  the  rails.  Intuitively,  we  expect  there  to  be 
concentrations  of  magnetic  field  at  the  rear  corners  of  the  interface  liquid,  with 
the  force  reaching  relatively  high  values,  perhaps  sufficient  to  eject  the  molten 
layer  forward  of  the  accelerating  mass  of  the  armature.  The  main  purpose  of  this 
report  is  to  show  that  this  is  not  the  case.  Let  dF{x)/dx  denote  the  gradient  of 
the  average  EM  force  over  the  armature  surface,  where  x  is  distance  measured 
between  the  rails.  We  shall  use  the  Biot-Savart-Lorentz  (BSL)  formula  to  deduce 
an  analytic  expression  for  this  derivative  and  from  this,  show  that  there  is  only 
a  small  variation  in  the  electromagnetic  force  over  the  rear  of  the  armature,  30% 
being  typical. 

Integrating  our  expression  for  the  force  gradient  to  obtain  the  total  force  F 
acting  on  the  armature,  and  using  the  standard  expression,  F  =  \L’F  where 
I  is  the  total  electric  current  and  L'  is  the  gradient  of  the  self-inductance,  we 
deduce  an  analytic  expression  for  L'  as  a  function  of  =  Kfwr  and  s  =  tfwr, 
where  br  is  the  distance  between  rails  of  width  Wr  and  thickness  t.  This  is  not 
exact,  since  we  have  assumed  the  current  distribution  to  be  uniform.  But  for 
small  values  of  5,  this  appears  to  be  a  good  approximation.  Our  value  for  L'  at 
^  —  0.78,  s  =  0.2  (the  lAT  rail  gun  ratios)  is  0.486  /iH/m,  in  good  agreement 
with  observation.  Using  numerical  methods,  Kerrisk  [1]  obtained  a  table  of  values 
that  for  the  given  parameters,  yields  L'  =  0.457  /iH/m. 

Finally  we  investigate  the  claim  that  the  original  Amperian  law  for  the  force 
between  two  current  elements  (which  is  directed  along  the  vector  joining  the 
elements),  gives  ‘very  nearly  the  same  acceleration  force’  as  the  BSL  law.  We  find 
that,  typically,  the  Amperian  law  overestimates  this  force  by  50%,  so  that,  apart 
from  lacking  a  physical  mechanism,  this ‘law  has  no  merit  from  an  experimental 
point  of  view.  It’s  prediction  about  the  buckling  of  rails  by  the  recoil  forces  is 
wrong. 

1.  The  Biot-Savart-Lorentz  Force  La'w 

The  force  element  dFi  acting  on  a  current  /j  dsi  due  to  a  current  I2ds2  at  a 
distance  r  (e.g.  see  [2])  is  given  by 

=  ^^^{‘^SiX(ds2Xr)}.  (1) 

This  relation  is  sometimes  referred  to  as  Ampere’s  law,  although  his  version  made 
dFi  parallel  to  r  (see  §5). 

The  configuration  shown  in  Fig.  1,  in  which  a  sheet  current  passes  down  each 
side  of  rail  2,  then  across  the  armature  surface,  and  up  each  side  of  rail  1,  shows 
our  simplification  of  the  geometry.  The  current  flowing  on  the  edges  of  the  rails 
is  ignoied,  oi  rather  divided  equally  between  the  inside  and  outside  surfaces.  To 
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Figure  1:  Sheet  Current  down  Rail  2,  across  the  Armature  and  up  Rail  2 


simplify  the  presentation,  we  shall  start  by  supposing  the  rails  to  have  vanishing 
thickness,  and  then  later  modify  our  solution  to  allow  for  the  rail  thickness,  t. 
Our  main  assumption  is  that  the  current  is  uniform,  which  is  reasonably  accurate 
for  small  values  of  the  ratio  s,  defined  above. 

The  armature  surface  carrying  the  current  between  the  rails  lies  on  y  =  0  in 
—a  <  X  <  a,  —b  <  z  <  b.  The  rails  lie  on  x  =  ±a,  0  <  y  <  oo,  —b<z<b.  We 
shall  consider  the  force  acting  on  an  element  of  current  {I /2b)  dzi  dx.  at  (x,  0,  zi) 
on  the  armature,  due  to  current  elements  (I /2b)  dz2dy  at  (±a,  y,  Z2).  Adopting 
unit  vectors  x,  y,  z,  we  have  for  rail  1  (see  Fig.  2), 


dx  =  x(ix,  dy  =  ydy,  r  = -(a  -  x)x  -  yy  +  (zi  -  Z2)z, 


hence  dxx(dyxr)  =  —  (a  — x)  dxdyy.  A  similar  expression  holds  for  rail  2.  Then 
(1)  yields  the  element  of  force 


dF{x,y,zi,Z2)  = 


in  the  y-direction. 


j _ a-  X _ 

1  [(a  —  x)2  +  y^  +  (zi  —  ^2)^]^ 


+ 


a  +  X 

[(a  +  x)2  +  y2  +  (zi  -  Z2)2]^ 


I 


dx  dy  dzi  dz2  , 


(2) 


First  we  shall  remove  the  Z\  and  Z2  dependence  by  integrating  over  —b<zi<b 
and  —b<Z2<  b.  Let  gi  =  (a  —  x),  §2  =  (a  -1-  x)  and  first  integrate  over 
—b  <  Zi  <  b.  It  follows  from  (2)  that  we  need 


/6-Z2 
-6— Z2 


9i  dC 


-b-z,  [y2  +  y2  + 

9l{i> 


>-) 


(C  =  Zl  -  22,  rfC  =  dZi). 

9i  {b  +  Z2) 


+ 


{9i  +  ?/^)bi  +  +  {b-  22)^]^  (5?  +  +y'^  +  {b  +  22)2] ^ 


2) 


Next  we  integrate  over  -b  <  Z2  <  b.  The  first  term  in  the  above  expression  yields 

and  a  similar  formula  results  from  the  second  term.  In  this  way  equation  (2)  is 
integrated  over  Zi  and  Z2  to  yield 

U  P 

dF{x,y)  =  -^{Fi(a:,y;^i)  +  F2(x,y;p2)},  (3) 


where 


Fi{x,y,gi) 


(pi  +y‘^  +  46^)^ 


(Pi +y^)^ 


and  similarly  for  F2(x,  y;  ^2)- 


2.  The  Force  Distribution  over  the  Armature 

The  next  integration  is  over  0  <  y  <  00,  where  we  are  assuming  that  the  rails 
are  long  compared  with  the  dimensions  a  and  b  shown  in  the  figures.  We  could,  at 
the  cost  of  some  algebra,  remove  this  restriction,  and  obtain  expressions  for  the 
force  on  the  armature  at  the  earliest  stages  of  the  acceleration  of  the  armature. 
This  might  be  worth  the  effort  at  a  later  stage. 

To  integrate  Fi{x,  y),  we  need  to  evaluate  integrals  of  the  type 


Figure  2;  Interacting  Current  Elements 
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Figure  3:  Distribution  of  the  Force  Ratio 


The  first  integral  is  a  standard  logarithmic  type,  while  the  second  is  transformed 
into  a  ‘sin“^’  type  by  the  substitution  2:  =  l/(j/^  +  The  outcome  is 

00  2  "  46^ 

r dF,{x, V, gi)dy^  I51  + 1-  i’f  +  =‘”'*(462  +W 


10  ~  '^i+462> 

Thus  we  arrive  at  the  required  distribution 


dF{x)  = 


/io/' 


SttP 
+Ua  +  x)  In 


fcTT  +  5  (a  —  x)  In 


(a  —  xy 


2,^  (a  _  3.)2  4.  4^2 

(a  +  xY 


+  6sin“ 


452  +  ^2. 

1  46^  —  (a  —  xY 
4b^  +  (a  —  xY 


_i  46^  -  (a  +  xY 

(a  +  xY  +  462  +  462  +  (a  +  x)2j  ’ 


(4) 


where  we  have  reversed  the  sign  to  obtain  the  force  magnitude  Expressed  in  a 
non-dimensional  form,  this  distribution  reads 


dF{x/a)  = 


TT  -I-  rln 


1  +  r' 


+  sin  ^  ^  -f-  s  In  ^ 


1  +  r^ 


1  +  s^ 


+  sin 


-1 


1  — 
1  -H 


(5) 


where 


_  1  —  x/a 
""  26/a 


_  1  -f  x/a 
^  ~  26/a 


We  have  evaluated  the  ratio  R  =  dF{x/a)/dF{0)  in  the  range  0  <  x  <  a. 
This  is  shown  in  Fig.  3,  for  the  choice  it  =  a/6  =  K/wr  =  0.78.  The  ratio  climbs 
from  unity  at  x  =  0  to  1.3058  at  the  edge  of  the  armature,  x  =  a,  a  surprisingly 
modest  increase. 

3.  The  Inductance  Gradient 

To  obtain  the  total  force  acting  on  the  armature,  we  integrate  (4)  over  —a  < 
X  <  a.  We  shall  allow  for  the  thickness  t  of  the  rail  at  a  later  stage,  i.e.  referring 
to  Fig.  1,  we  shall  start  with  the  case  a  =  ai,  as  depicted  in  Fig.  2.  Two  types  of 
integral  appear.  The  first  is 
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and  the  second  is 


f 

J  —a 


2  _  1 

dx  =  -267^sin~^^^^  ^  +  261n(7j^  +  1) , 


46^  +  gi 

where  we  have  introduced  the  ratio  n  =  afb.  Thus,  evaluating 

/  dF{x)  dx  =  F , 

J —a 

we  arrive  at  the  EM  force  acting  on  the  armature: 

F  =  -  -  ln(^^)  +  i  ln(7t2  +  1)  -  ?  sin-‘ 

-1  I  TT^-R^^Tm  -K  + 


(7) 


(8) 


In  the  above  model,  we  have  placed  all  the  current  in  the  rail  on  a  single  sheet. 
Our  final  calculation  is  to  take  the  thickness  of  the  rails  into  account.  The  same 
type  of  integrals  appear  as  above,  except  that  the  contribution  from  the  current 
sheet  on  the  outside  of  the  rails  requires  us  to  replace  a  by  Oi  (see  Fig.  1).  Thus 
for  this  sheet,  instead  of  (6),  we  have  an  integral  of  the  same  form,  except  that, 
whereas  in  (6),  =  a  -  x,  we  now  have  gi=ai-  x.  As  the  range  of  the  integral 

remains  -a  <  x  <  a,  this  change  merely  complicates  the  algebra.  We  find  that 


/  1^1  In 


9i 


gi  +  4b^ 


dx  =  I  (7^^  +  5^)  ln( 


+  1  + 
V?  +  5^ 


) 


+  1+^2  1  +s2 

+  ln( 


where  5  =  (oi  -  a)/2b.  Similarly  (7)  is  generalised  to 
r  sin-‘  dx  =  -26  ((,;  +  .) 

J-a  462  +  p2  I^V  !  V(^  +  5)2  +  1; 


(9) 


I  _  ^2  \ 

+  +  .  (10) 

There  are  also  the  integrals  for  the  opposite  side  of  the  rail  system,  i.e.  where 
52  =  02  +  X. 

We  shall  add  the  contributions  to  F  from  the  inside  and  outside  of  the  rails, 
using  weighting  factors  1  -  and  ^f.  Thus  at  /  =  1,  each  surface  contributes  the 
amount  to  F  that  a  uniform  current  distribution  would  require.  In  the  absence 
of  details  about  the  distribution,  this  is  the  least  biased  choice.  We  use  the 
relation  F  =  to  deduced  the  value  of  U .  Collecting  the  results  together, 

and  adopting  the  dimensionless  numbers 


R'  —  bj-jiVj’  s  —  i/jiVj- , 


(11) 
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Figure  4:  Rail  dimensions 


where  br  is  the  distance  between  rails  of  width  Wr  and  thickness  t  (see  Fig.  4), 
our  final  expression  for  L',  in  fxH/m,  is 


L'  =  y%{\  -  (2  -  /)4i  -  /O.},  (0  <  /  <  1) , 


(12) 


where 


and 


^  1  f(7J  +  5)^,  /(rc4-5)^  +  l\  1  /(■R  +  5)^  +  l\  2  \2  ,  1  1 

$0  =  -{- - ^ln(^ - +  -ln(  ,,  ,  xo--)  --ln{(7i  +  5)"  +  l} 

27r  \  7^  (tc  +  5)2  /  TC  (1  +  5)2  >'  TZ 


+2- 


(7^  +  5)^^^_l/(7^  +  5)2-l^  25  ,/l^ 


7^ 


Sin 


( 


(7^4' 5)2  +  1^  n 


5^  /l+5^\  2  .  2\ 


Table  1:  Values  of  L'{n,s) 


5\^ 

0.4 

0.6 

0.8 

1.0 

1.5 

2.0 

0 

0.395 

0.502 

0.577 

0.628 

0.692 

0.699 

0.01 

0.344 

0.457 

0.548 

0.624 

0.771 

0.880 

0.10 

0.318 

0.426 

0.513 

0.587 

0.729 

0.836 

0.15 

0.305 

0.410 

0.496 

0.568 

0.708 

0.814 

0.20 

0.293 

0.396 

0.480 

0.550 

0.689 

0.793 

0.25 

0.282 

0.382 

0.465 

0.534 

0.671 

0.774 

0.30 

0.271 

0.370 

0.451 

0.519 

0.654 

0756 

0.50 

0.236 

0.328 

0.404 

0.468 

0.597 

0.695 
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The  above  table  sets  out  values  of  the  inductance  gradient  as  functions  of  n 
and  s.  The  following  table  sets  out  some  values  taken  from  Kerrisk’s  paper  [1]. 
The  bracketed  values  are  deduced  from  (12).  The  agreement  is  fair  for  s  =  0.1 
(thin  rails).  It  remains  uncertain  whether  the  discrepancies  at  larger  values  of 
s  are  due  to  the  approximations  introduced  by  Kerrisk  or  to  the  assumption  of 
uniform  current  density  adopted  in  this  paper. 


Table  2:  Kerrisk’s  values  for  V 


0.2 

0.5 

1.0 

1.5 

3.0 

0.10 

0.183 

0.361 

0.555 

0.690 

0.942 

(0.182) 

(0.375) 

(0.587) 

(0.729) 

(0.991) 

0.50 

0.171 

0.325 

0.491 

0.607 

0.832 

(0.121) 

(0.284) 

(0.468) 

(0.597) 

(0.840) 

4.  Ampere’s  law, 

Ampere’s  original  formula  was  based  on  the  assumption  that  the  force  between 
two  current  elements  was  directed  along  the  line  joining  them.  The  main  weakness 
of  his  plausible  treatment  was  that  no  mechanism  for  the  transfer  offeree  between 
the  elements  was  advanced,  which  gave  Maxwell  some  concern,  although  at  the 
end  of  his  own  investigations,  he  was  content  to  accept  Ampere’s  description. 

In  a  letter  to  William  Thomson  dated  13th  Nov.,  1854  [3]  he  wrote; 

I  got  up  the  fundamental  principles  of  electricity  of  tension  easily 
enough.  I  was  greatly  aided  by  the  analogy  of  the  conduction  of  heat, 
wh:  I  believe  is  your  invention  at  least  I  have  never  found  it  elsewhere. 

But  then  I  tried  to  make  out  the  theory  of  the  attractions  of  currents 
but  tho  I  could  see  how  the  effects  could  be  determined  I  was  not  sat¬ 
isfied  with  the  form  of  the  theory  which  treats  of  elementary  currents 
&:  their  reciprocal  actions,  &  I  did  not  see  how  any  general  theory 
was  to  be  formed  from  it.  I  read  Ampere’s  investigations  this  term 
&  greatly  admired  them  but  thought  there  was  a  kind  of  ostensive 
demonstration  about  them  wh:  must  have  been  got  up,  after  Ampere 
had  convinced  himself,  in  order  to  suit  his  views  of  philosophical  in¬ 
quiry,  and  as  an  example  of  what  it  ought  to  be.  Yet  I  believe  that 
there  is  no  doubt  that  Ampere  discovered  the  laws  k  probably  by  the 
method  wh:  he  has  given. 

In  a  letter  dated  28th  Dec.,  1861,  written  to  Henry  Richmond  Droop,  he 
returns  to  the  subject  with: 

I  am  trying  to  form  an  exact  mathematical  expression  for  all  that  is 
known  about  electro-magnetism  without  the  aid  of  hypothesis,  and 
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also  what  variations  of  Ampere’s  formula  are  possible  without  con¬ 
tradicting  his  expressions.  All  that  we  know  is  about  the  action  of 
closed  currents — that  is  currents  through  closed  curves.  Now  if  you 
make  a  hypothesis  (1)  about  the  mutual  action  of  the  elements  of  two 
currents,  and  find  it  agree  with  experiment  on  closed  circuits,  it  is  not 
proven,  for — 

For  if  you  make  another  hypothesis  (2)  which  would  give  no  action 
between  an  element  and  a  closed  circuit,  you  may  make  a  combination 
of  (1)  and  (2)  which  would  give  the  same  result  as  (1).  So  I  am 
investigating  the  most  general  hypothesis  about  the  mutual  action 
of  elements,  which  fulfils  the  condition  that  the  action  between  an 
element  and  a  closed  circuit  is  null. 

Maxwell  wrote  to  Droop  again  on  24th  Jan.,  1862: 

I  want  to  see  if  there  is  any  evidence  from  the  mathematical  expres¬ 
sions  as  to  whether  element  acts  on  element,  or  whether  a  current 
first  produces  a  certain  effect  in  the  surrounding  field,  which  after¬ 
wards  acts  on  any  other  current.  ... 

As  a  fact,  the  effect  on  a  current  at  a  given  place  depends  soley  on  the 
direction  and  magnitude  of  the  magnetic  force  at  that  point,  whether 
the  magnetic  force  arises  from  currents  or  from  magnets.  So  that  the 
theory  of  the  effect  taking  place  through  the  intervention  of  a  medium 
is  consistent  with  fact,  and  (to  me)  appears  the  simplest  in  expression 


It  appears  that  Maxwell  was  very  close  to  deducing  the  law  that  much  later 
was  to  be  attributed  to  Lorentz.  His  uncertainty  about  whether  or  not  it  was 
valid  to  write  an  expression  for  the  force  between  current  elements  is  evident 
from  the  above  quotations. 

The  consensus  of  scientific  opinion  today  is  that  Ampere’s  original  form  is 
not  correct,  although  Graneau  [4,5,6,7]  has  persisted  with  the  Amperian  model, 
and  has  collected  various  experimental  results  that  appear  to  support  his  view. 
It  appears  that  his  principal  aim  has  been  to  locate  the  recoil  force  in  the  rails 
close  to  the  armature.  Figure  5  illustrates  the  difference  between  the  two  laws 
and  makes  clear  that,  if  he  is  correct,  there  would  be  very  large  buckling  forces 
acting  on  the  rails  just  behind  the  accelerating  armature.  In  his  1982  paper  he 
concluded: 

Longitudinal  Ampere  forces  reveal  not  only  an  unexpected  recoil 
mechanism,  but  they  also  suggest  that  the  projectile  branch  is  be¬ 
ing  stressed  like  a  strut  in  compression.  The  recoil  forces  should  be 
capable  of  bucking  copper  rails  unless  the  rails  are  suitably  reinforced. 
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The  broad  conclusion  drawn  from  this  investigation  is  that  Ampere’s 
law  correctly  describes  the  distribution  of  mechanical  forces  of  electro¬ 
magnetic  origin  around  the  railgun  circuit.  It  supports  the  intuitive 
notion  of  the  recoil  force  residing  in  the  rails  right  adjacent  to  the 
accelerating  projectile.  Both  laws  [Lorentz  plus  Ampere]  give  very 
nearly  the  same  acceleration  force,  and  they  also  agree  on  the  distri¬ 
bution  of  this  force.  But  the  Amperian  law  also  predicts  an  additional 
compressive  force  to  be  found  in  the  projectile  branch  which  does  not 
conform  with  the  Lorentz  force  concept. 
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(b)  Force  Distribution  According  to  Lorentz's  Law 

Figure  5:  Force  Distribution  due  to  (a)  Ampere’s  law  and  (b)  Lorentz’s  law. 

In  his  1987  articles  Graneau  cites  an  experiment  showing  the  buckling  of  thin 
aluminium  rails.  While  there  can  be  no  doubt  about  the  observations,  it  does 
not  follow  that  the  Amperian  force  must  necessarily  exist— there  may  be  other 
explanations  of  the  buckling.  There  are  two  problems  with  Graneau’s  appeal 
to  a  discarded  nineteenth  century  law;  first  there  is  the  question  of  macroscopic 
mechanism  raised  by  Maxwell,  (see  his  second  letter  to  Droop  quoted  above) 
and  secondly  no  microscopic  mechanism  has  been  described,  i.e.  how  do  the  elec¬ 
trons  in  the  rail  near  the  armature  respond  to  the  motion  of  the  electrons  in  the 
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In  his  seminal  treatise  on  Electricity  and  Magnetism  [8],  Maxwell  made  a 
complete  analysis  of  the  possible  force  laws  between  current  elements.  He  arrived 
at  four  possible  choices,  each  one  of  which  was  consistent  with  the  closed  circuit 
constraints.  On  page  161  he  writes 

Thus  Ampere  assumes  that  the  force  between  two  elements  is  in  the 
line  joining  them.  This  gives... 

Grassmann  assumes  that  two  elements  in  the  same  straight  line  have 
no  mutual  action.  This  gives... 


In  fact  Grassmann’s  empirical  choice  is  equivalent  to  the  BSL  law  in  equation 
(1).  After  dealing  with  the  two  remaining  and  unlikely  possibilities,  Maxwell 
concludes: 

Of  these  four  different  assumptions  that  of  Ampere  is  undoubtably 
the  best,  since  it  is  the  only  one  that  makes  the  forces  on  the  two 
elements  not  only  equal  and  opposite  but  in  the  straight  line  which 
joins  them. 

The  greatest  19th  Century  scientist  made  very  few  errors,  but  this  was  one  of 
them. 

The  assertion  that  both  laws  give  very  nearly  the  same  acceleration  force  is 
one  that  we  can  readily  check  with  some  analysis  similar  to  that  given  above  for 
the  BSL  law. 


5.  Inductance  Gradient  According  to  Ampere’s  Law, 

Maxwell  shows  that  the  force  element  (on  either  current) 

fJ'ohh 


dF 


or 


dF 


iJ-Jxh 


(3f  •  dsi  f  •  ds2  +  2dsi  •  (is2)f , 
(3  cos  6i  cos  92  +  2  cos  e)  r, 


(13) 


where  the  angles  9i,  62  and  e  are  shown  in  Fig.  6.  Referring  to  Fig.  2,  we  find 
cos  01  =  {x  —  a)/r,  cos  02  =  J/Ai  ^  =  '^/2,  and  hence  in  place  of  (2)  we  get 


dF{x,y,zi,Z2)  = 


3/io/2 


y{a-  x) 


\  [(a  -  x)2  +  t/2  +  (zi  -  22)^]^ 
the  omitted  term  having  (a  +  x)  in  place  of  (a  —  x). 

Integrating  over  0  <  y  <  00,  —a  <  x  <  a,  we  get 


+ 


dF{zuZ2)  = 


3noP  /4a^  +  (zi  -  Z2y 

3^  V  {^1-^2)^ 


dzx  dz2. 
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Figure  6:  Geometry  for  the  Amperian  Law 


The  final  integrations  are  straightforward.  Dividing  the  result  by  we  arrive 
at  the  inductance  gradient  (in  //H/m) 

—  7^^  ln(l  +  +  An  tan~^  1  /n),  (14) 

which  is  to  be  compared  with  (12)  at  /  =  0.  (We  shall  not  trouble  to  allow  for 
the  rail  thickness,  since  the  conclusion  we  wish  to  reach  does  not  require  this.) 


Table  3:  Comparison  of  Ampere’s  and  Lorentz’s  Inductance  Gradients 


7^ 

0.4 

0.6 

0.8 

10 

1.5 

2.0 

L'(A) 

0.521 

0.691 

0.828 

0.942 

1.164 

1.328 

L'(BSL) 

0.395 

0.502 

0.577 

0.628 

0.692 

0.699 

ratio 

1.32 

1.38 

1.44 

1.50 

1.68 

1.90 

It  will  be  seen  from  the  above  table  that  the  use  of  the  original  Amperian 
law  gives  values  for  the  inductance  gradient  very  much  larger  than  the  correct 
formula,  labelled  BSL.  And  since  experimental  measurements  of  the  acceleration 
yield  values  of  L'  fairly  close  to  the  BSL  values,  we  must  conclude  that  Graneau’s 
assertion  that  both  laws  give  very  nearly  the  same  acceleration  force  is  not  sup¬ 
ported.  It  also  follows  that  the  idea  that  the  recoil  forces  are  responsible  for 
buckling  of  rails  is  wrong. 
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Summary 


This  report  returns  to  the  melt-wave  model  described  in  earlier  reports  and 
is  concerned  with  the  question  ‘what  observations  could  be  made  to  test  it?’ 
Measurements  could  be  made  of  the  thickness  6m  of  the  aluminium  layer  deposited 
on  the  copper  rails  and  (with  some  difficulty)  of  the  depth  h  of  the  aluminium 
layer  removed  by  melting  from  the  armature.  Direct  measurement  of  the  melt- 
wave  speed  would  be  difficult,  so  we  are  left  with  6m  as  providing  the  only  link 
between  theory  and  experiment. 

In  [1],  the  depth  h  of  the  material  removed  from  the  armature  by  melting 
was  determined  from  the  diffusivity  of  the  electric  current  on  the  assumption 
that  the  mechanism  was  essentially  adiabatic.  In  [2]  doubts  about  the  validity  of 
this  were  advanced,  it  being  argued  that  the  thermal  diffusivity  should  have  the 
dominant  role.  During  the  author’s  visit  to  lAT  in  October,  1996,  the  variation 
of  6m  with  distance  x  along  the  rail  was  obtained  for  one  experiment.  This 
relationship  follows  different  laws  for  thermal  and  ohmic  diffusivity;  the  lAT 
experiment  clearly  favoured  thermal  diffusivity.  In  this  report  we  present  a  semi- 
empirical  theory  for  the  dependence  of  h  on  the  thermal  diffusivity,  the  time  of 
flight  of  the  armature  and  its  velocity.  From  this  an  expression  is  obtained  for 
the  thickness  6m,  including  the  effects  of  viscous  heating,  taken  from  [3].  The 
theoretical  values  this  thickness  are  found  to  be  in  good  agreement  with  the 
measured  values  for  the  experiment  in  question.  It  remains  to  test  the  theory 
with  other  observations. 

1.  Introduction 

In  the  melt-wave  model  [1]  there  are  four  variables  that  could  possibly  be 
measured  to  check  the  theory.  These  are  the  thickness  h  of  the  material  removed 
from  the  armature  by  the  wave,  the  thickness  6m  of  the  material  deposited  on  the 
copper  rail,  the  speed  of  the  melt  wave,  Vm  and  perhaps  the  time  t*  for  the  onset 
of  transition,  which  is  assumed  to  occur  when  the  wave  has  removed  sufficient 
material  to  interrupt  the  electrical  contact  between  the  armature  and  rail.  The 
first  three  of  these  variables  are  illustrated  in  Fig.  1.  The  difficulty  with  the  last 
variable  is  that  it  depends  on  the  force  between  the  rail  and  armature,  so  goes 
beyond  the  theory  of  the  melt  wave.  Also  the  melt-wave  speed  Vm  would  be 
difficult  to  measure,  which  leaves  us  with  h  and  6m-  When  viscosity  is  ignored, 
these  variables  are  related  by  continuity, 

PsVmh  =  peiV  +  Vm)6m  ,  (l) 

where  and  p(  are  the  densities  of  the  aluminium  in  the  solid  and  liquid  states, 
and  V  is  the  relative  speed  of  the  rail  and  armature. 

An  expression  for  Vm  can  be  deduced  from  the  energy  equation  applied  to  the 
solid  legion  of  the  armatuie  extending  in  front  of  the  melt  wave.  In  [1]  we  showed 
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Figure  1:  Melt- wave  model 


that 

^  _  2fx,GV 

2r]rPQw^  TTTJr 

where 

Q  =  Cj,{T^-To)  +  Lm,  G=^^%f^Xa[l  +  (— ft.  (3) 

W  'Xr 


In  these  expressions  Tm  is  the  melting  temperature,  To  is  the  laboratory  tem¬ 
perature,  T7„  is  resistivity  in  the  solid  armature,  averaged  over  {Tm,  To),  i]r  is  the 
rail  resistivity,  /  is  the  total  current  flowing  through  the  armature,  /Xq  is  the  per¬ 
meability  of  free  space,  w  is  the  depth  of  the  surface  of  the  armature  in  contact 
with  a  rail,  Lm  is  the  enthalpy  of  fusion  (latent  heat),  Cp  is  the  specific  heat  at 
constant  pressure,  and  Xa,  Xr  are  the  thermal  diffusivities  in  the  armature  and 
rail.  Note  that  Vm  is  independent  of  the  size  of  the  gap  thickness,  h. 

In  [1]  we  also  made  a  case  for  the  choice 

TTTJr 
2^loV’ 

this  being  the  length  scale  (in  the  direction  of  motion)  for  the  diffusion  of  the 
electric  current,  concentrated  by  the  relative  motion  of  the  rail  and  armature. 
With  this  model,  thermal  diffusivity  plays  no  role  in  determining  h.  It  was  based 
on  the  notion  that  the  mechanism  is  essentially  adiabatic,  i.e.  that  the  thickness 
of  the  gap  is  not  the  (negligible)  length  scale  for  the  penetration  of  the  heat,  but 
the  distance  over  which  melting  occurs,  and  this,  it  was  supposed,  is  limited  to  the 
region  of  the  most  intense  heating — ohmic  heating.  However  in  [2]  some  doubts 
were  expressed.  And  during  the  author’s  last  visit  to  lAT  (Sept.,  1996),  the 
experimental  evidence  supported  a  thermal  rather  than  an  electric  explanation  of 
h.  In  [3]  we  considered  the  effect  of  viscous  heating,  which,  being  dispersed  over 
the  length  of  the  armature,  has  little  influence  on  h. 


(4) 
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In  this  report  we  shall  make  three  changes  to  the  theory:  first  we  shall  give  a 
thermal  treatment  of  the  gap  thickness  h,  secondly  we  shall  reconsider  the  loses 
due  to  thermal  conductivity,  represented  by  the  parameter  G  in  (3)  and  thirdly 
we  shall  include  the  effect  of  viscous  heating  on  6m  ■  This  will  allow  us  to  obtain 
expression  for  6m  that  is  much  more  accurate  than  that  obtained  from  (1)  to  (4). 
Finally  our  formula  will  be  compared  with  some  observations.  This  appears  to 
be  the  only  test  we  can  make  of  the  melt-wave  model  at  this  stage. 

2.  The  Spread  of  Current  in  the  Armature 

The  current  at  the  interface  is  concentrated  into  a  strip  only  6r  =  rir/l^aV  « 
2.12xl0“^/y  wide,  for  example,  at  a  value  of  1^  =  330m/s,  6r  —  64  ^m.  With  the 
distance  x  measured  from  the  edge  of  the  melt  wave,  in  the  direction  of  motion 
of  the  armature,  in  [1]  it  is  shown  that  the  current  density  across  the  interface 
has  the  the  distribution 

jy  =  (-00  <  X  <  0),  (jo  =  const.) .  (5) 

The  physical  argument  for  the  length  scale  6r  is  as  follows.  In  a  time  t  the 
field  in  the  rail  diffuses  through  a  distance  {^rt)K  where  (r  =  ?7r//^o  is  the  electric 
diffusivity.  In  the  same  time,  the  motion  of  the  rail  relative  to  the  armature,  along 
the  positive  OX-axis,  carries  the  rail  a  distance  Vt,  and  if  these  displacements 
are  opposite  and  cancelling  so  as  to  produce  a  steady  state,  we  can  eliminate  t 
to  obtain  6r  —  G/V ,  as  required.  This  is  what  happens  in  the  rail  adjacent  to 
the  armature,  viewed  as  being  stationary— as  fast  as  the  field  diffuses  forwards, 
it  is  swept  back  in  the  rail.  The  boundary  condition  at  the  rail/armature  inter¬ 
face  imposes  the  same  OX-scale  length  on  the  current  flowing  normally  into  the 
armature,  i.e.  jy  varies  in  the  OX-direction  on  a  length  scale  of  6r  in  both  the 
rail  and  the  armature. 

However  once  the  current  enters  the  armature  and  is  not  longer  being  swept 
along  OX  by  the  motion  of  the  rail,  it  spreads  out  laterally,  and  in  t  seconds 
reaches  a  distance  ~  0.292  fi,  where  we  have  used  the  armature 

value  =  10.7x10  ®  (see  below).  Thus  in  just  10 /zs,  Aq  ~  1mm,  which  is 
about  fifteen  times  the  scale  length  6r  at  the  interface.  We  have  depicted  this 
phenomenon  in  Fig.  2. 

Notice  the  concentration  of  current  at  the  interface  just  in  front  of  the  gap;  it  is 
rather  like  a  point  source,  from  which  heat  will  spread  radially  into  the  armature. 
The  rail  remains  relatively  cool  because  in  a  frame  fixed  in  the  armature,  it  is 
moving  rapidly.  We  shall  return  to  the  question  of  the  heat  flow  into  the  rail 
shortly.  We  have  marked  a  circular  region  of  radius  h,  centred  on  the  ‘heat 
souice  ,  S,  and  extending  in  the  direction  of  motion  and  into  the  armature.  We 
shall  jissume  that  the  armature  material  in  this  region  has  just  been  melted  and 
is  on  the  point  of  being  entrained  on  the  rail  and  thus  removed  from  the  region 
of  interest.  We  shall  deal  with  the  problem  of  determining  h,  after  considering 
the  tluirmal  losses  into  the  rail. 
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Figure  2:  Spread  of  Current  in  the  Armature 


3.  Thermal  Losses 

The  scale  length  for  heat  diffusing  from  a  source  into  a  stationary  medium  of 
thermal  diffusivity  x,  is  A  =  {xt)  ^ ,  where  t  is  the  time  from  the  moment  the  heat 
from  the  source  is  liberated.  With  both  copper  and  aluminium,  x  ^ 

Thus  after,  a  typical  time  lapse,  when  the  armature  is  about  half  way  towards 
the  muzzle,  say  at  t  =  1ms,  A  is  about  316 /im.  This  means  that  at  distances 
greater  than  this  from  the  current  front,  the  temperature  has  not  been  affected 
by  ohmic  heating. 

In  [1]  we  assumed  that  the  thermal  loss  term  defined  as  G  in  (3),  viz. 


was  a  satisfactory  approximation.  However  the  term  (Xo/Xr)^  needs  to  be  re¬ 
considered.  It  appeared  as  the  ratio  Aa/Ar,  where,  for  a  given  time  interval  t, 
Xa  is  the  thermal  skin  depth  (XoO^  iii  th®  armature  and  Ar  is  the  thermal  skin 
depth,  wrongly  assumed  to  be  (Xrt)^  in  the  rail.  The  heat  losses  were  taken  to 
be  transverse  only,  i.e.  in  the  OF-direction  into  the  rail  and  armature.  There  is 
very  little  difference  between  Xa  and  Xr,  so  we  set  the  ratio  Ao/Ar  equal  to  unity. 
This  is  wrong.  The  failing  here  is  that  the  motion  of  the  rail  has  been  ignored. 
That  heat  flux  into  the  rail  is  affected  by  the  relative  motion  of  the  heat  source 
and  the  rail  being  heated  was  overlooked. 

The  classical  treatment  of  the  problem  just  described  is  to  be  found  in  [4]; 
also  there  are  some  interesting  applications  to  welding  and  similar  problems,  with 
moving  heat  sources,  in  [5].  The  situation  is  analogous  to  that  described  in  §2  for 
the  physical  origin  of  6^.  In  the  direction  of  motion,  i.e.  along  —OX,  the  diffusion 
length  in  the  rail  is  Xr/V-  Since  Xr  ~  Xa)  we  shall  ignore  the  small  difference 
between  these  diffusivities  and  write  them  both  as  x- 
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Figure  3:  Diffusion  of  Heat  into  the  Rail 


Consider  the  trajectory  of  a  heat  spot,  H,  moving  into  the  rail,  orthogonally 
across  the  interface  (see  Fig.  3).  In  a  time  t,  H  will  move  in  the  OY-direction 
a  distance  y  =  and  in  the  same  time,  it  will  be  displaced  along  OX  a 

distance  Vt.  Eliminating  t,  we  obtain  the  trajectory  =  {x/V){xo  -  x),  where 
xq  is  the  point  on  the  interface  at  which  the  spot  enters  the  rail.  Referring  to  Fig. 
3,  we  see  that  =  xq  =  {xt'f  and  =  y  at  a;  =  0,  i.e.  A,.  = 

We  therefore  obtain  the  ratio 


^ 

Ar  “  x^V^ 


(6) 


For  example,  at  t  =  1  ms,  1/  =  330  m/s,  X^/X^  «  32.  Previously,  we  set  this  ratio 
equal  to  unity,  and  thus  considerably  underestimating  the  heat  losses. 


4.  The  Gap  Thickness 

It  was  assumed  in  [1]  that  the  melt  wave  had  a  depth  h  equal  to  7r^,./2. 
This  was  on  the  grounds  that  (i)  the  melting  occurred  where  the  current  was 
most  concentrated  and  (ii)  that  the  gap  extended  to  the  point  where  the  normal 
current,  in  a  rectangular,  steady-state  model,  fell  to  zero.  The  first  point  is 
plausible,  except  that  over  the  range  — 7r5r/2  <  x  <  0,  the  current  density  in 
(5),  varies  from  a  maximum  value  of  jo  to  ~  0.21  joi  which  is  evidently  much  too 
large  a  range  to  encompass  only  that  part  of  the  armature  that  is  at  the  melting 
point.  The  second  point  is  not  valid,  since  the  current  distribution  obtained  in 
[1]  is  based  on  the  assumption  that  the  flow  is  in  a  steady  state  simultaneously 
in  both  the  rail  and  the  armature.  The  spread  of  current  into  the  armature,  as 
depicted  in  Fig.  2,  was  not  taken  into  account. 

The  problem  of  accurately  determining  the  size  of  the  melting  region  shaded 
in  Fig.  2  is  quite  complicated— for  an  accurate  determination,  numerical  methods 
are  unavoidable,  but  we  can  arrive  at  an  estimate  as  follows. 

First  we  note  that  heat  conduction  problems  involving  change  of  state  have 
been  extensively  studied,  especially  in  problems  involving  melting  and  freezing 
[4].  In  the  case  of  cylindrical  symmetry,  with  a  melted  region  r  <  R,  surrounded 
by  region  r  >  in  which  the  temperature  falls  to  a  value  Tq  at  infinity,  it  can  be 
shown  that  the  radius  R  of  the  region  at  the  melting  temperature  is  given  by 
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(see  p.  295  of  [4]) 


R  =  2i{xU)K 


where  i  is  the  root  of 


«V’Ei(-£^)  +  ^2i%-^  =  0,  (7) 

Lijn 

in  which  Ei(a:)  is  the  exponential  integral  and  the  time  is  the  time  measured 
from  the  instant  that  the  temperature  reaches  T,n  at  r  =  0.  With  aluminium  we 
find  that  i  =  0.494. 

However  the  value  0.988(xto)^  is  not  correct  for  the  radius  h  of  the  melted 
region  shown  in  Fig.  2,  since  it  does  not  allow  for  the  great  difference  between  the 
thermal  diflFusion  length  into  the  armature,  viz.  Aq  ~~  {xio)K  and  the  diffusion 
length  into  the  rail,  which  we  showed  in  §3  to  be  Ar  =  {(x/^)(x*o)^}^-  The  size 
of  the  melted  region  will  be  considerably  smaller  than  Aa,  since  at  least  half  of  the 
heat  from  the  ohmic  heat  source  will  not  be  available  for  melting  the  armature. 
In  fact  the  singularity  S  in  Fig.  2  is  more  like  a  dipole  than  a  heat  source,  with 
the  rail  providing  the  heat  sink.  The  effective  radius  h  of  the  melted  region  will 
lie  somewhere  between  the  values  Aa  and  Ar.  On  the  one  hand  its  value  must 
tend  to  Aa  as  V  tends  to  zero  and  on  the  other  hand,  it  must  vanish  as  the  rail 
velocity  tends  to  infinity,  since  in  this  case  the  heat  sink  is  dominant.  Notice 
that  Ar  is  the  geometric  mean  of  the  two  (orthogonal)  diffusive  lengths,  (x^o)^ 
and  x/V-  This  and  the  limiting  conditions  just  noted,  suggests  that  h  =  (ArAn)^ 
would  be  a  plausible  estimate  for  the  size  of  the  melted  region.  In  the  absence  of 
a  complete  theory,  we  shall  therefore  adopt  this  value. 

Hence,  with  the  value  x  ~  l  we  obtain 

h  =  {(x/lO^(x<o)*(x4^}^  =  (x^o)^(x/V")^  =  3.36x10-^4 /Vi  .  (8) 

We  regard  this  result  as  being  ‘semi-empirical’ — its  merit  to  be  checked  against 
experiment. 

5.  The  Thickness  of  the  Deposited  Layer 

The  first  step  is  to  modify  equation  (2)  by  replacing  (Xa/Xr)^  in  G  by  the 
ratio  given  in  (6) .  With  numerical  values  substituted,  we  get 

y  =  -  2.96xl0-^[l  -h  lOtivi] .  (9) 

However,  there  is  an  assumption  involving  the  value  of  the  average  fj^  that  needs 
to  be  revised.  In  the  derivation  of  (2),  we  wrote 

f  Vajo  dx  ft!  I  jl  dx  , 

•/—oo  J —oo 
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and  later  assumed  that  77^  was  the  arithmetical  average  of  the  limiting  values, 
lao  —  3,7x10“®  (25°C)  and  rjam  =  10.7xi0“®  (660°C).  But  the  correct  average  is 
evidently 

Va  = 

If  we  now  adopt  the  linear  approximation 

Va  =  Vao{l  +  0:{T  -  To)}  =  T]ao{l  +  a{Tm  -  , 

where  Aq  is  the  thermal  diffusion  distance,  we  get 

s  1 

Va  =  Tl—'naO  +  =  <^r/2Aa)  • 

1  +  s  1  +  s 


fO 

/  Vae^dy  {y  =  2x/6r). 
J  — CXD 


In  §2  we  obtained  the  typical  values  6^  =  64  /xm,  Aq  =  1000  jiva..  And  it  is 
easily  shown  that,  in  general,  s  <  1.  Therefore  we  may  write  v^  =  r/am-  This 
modification  replaces  \{r)aQ-^Vam)  by  Vam-  Hence  the  factor  1.04  in  (9)  is  replaced 
by  1.55.  With  this  change  and  the  values  w  =  18.85  mm,  /  =  10“®/,  i  =  10~H, 
(9)  becomes 

y  =  4.36xlO“®P  -  2.96xl0“®[l  +  ,  (10) 

where  the  current  is  in  kA  and  the  time  in  ms. 


From  (1) 


Hence  by  (8) 


dm  = 


Pe  1  +  Vm/V  1  +  Vm/V 


277  \ 


^/V 


l+Vm/V' 


where  to  is  now  in  ms  and  dm  is  in  /xm. 


(11) 


In  [3]  we  showed  that  viscosity  alone  contributes  a  thickness  0.59(ay)i  /xm 
to  the  deposited  layer.  Here  a  is  the  length  of  armature  in  contact  with  the 
rail,  measured  in  the  direction  of  motion.  But  with  the  melt  wave  continuously 
breaking  this  contact,  a  varies  between  the  total  length  a  and  a  much  shorter 
distance.  The  best  we  can  do  in  this  complicated  situation  is  to  put  a  =  la. 
With  0  =  1  cm,  the  viscous  thickness  is 


dmv  =  5.9xlO-Vi  .  (12) 

Fortunately  this  is  much  smaller  than  the  thickness  due  to  ohmic  dissipation  and 
therefore  its  errors  are  not  too  serious  in  the  total  thickness. 

The  total  thickness  is  therefore 


dm  =  277 


t!  Vm/V 


V'i  l+Vm/V 


+  5.9x10“  Vi  /xm. 


(13) 
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Table  1:  Comparison  of  Theoretical  and  Observed  values  of  Sm 


V  (m/s) 

341 

525 

647 

735 

800 

883 

941 

I  (kA) 

420 

410 

390 

377 

369 

356 

346 

t  (ms) 

0.471 

0.59 

0.68 

0.75 

0.81 

0.87 

0.93 

to  (ms) 

0.001 

0.12 

0.21 

0.28 

0.34 

0.40 

0.46 

VmlV 

0.69 

0.62 

0.54 

0.48 

0.45 

0.40 

0.36 

6mv  (Aim) 

1.09 

1.35 

1.50 

1.60 

1.67 

1.75 

1.81 

h  (/xm) 

4.8 

26 

31 

33 

35 

36 

37 

<5m  (mhi) 

3.1 

11.4 

12.2 

12.4 

12.4 

12.0 

11.7 

■5“"  im) 

3 

7 

10 

12 

11 

10 

10 

where  by  (10) 


Y  =  4.36x10"®/2  -  2.96xl0-^[l  +  .  (14) 

In  these  equations  to  and  t  are  in  ms  and  I  is  in  kA. 

The  distinction  between  t  and  to  is  important;  t  is  measured  from  the  begin¬ 
ning  of  the  armature’s  motion,  whereas  to  is  measured  from  the  instant  that  the 
armature  first  melts,  which  we  take  to  occur  when  6m  first  changes  from  zero. 

6.  Comparison  with  Observations 

To  test  the  theory,  we  shall  use  the  observations  of  lAT  experiment  95101202, 
in  which  an  aluminium  armature  of  mass  26.79  grms  was  accelerated  over  copper 
rails.  The  observed  values  appear  in  the  first  four  rows  of  the  following  table  and 
in  the  last  row.  Rows  five  to  eight  were  calculated  using  the  experimental  values 
of  the  first  four  rows  and  the  theory  given  in  §5. 

The  value  t  =  0.471  entered  in  the  first  column  was  chosen  to  give  the  ex¬ 
perimental  value  of  6m-  To  obtain  to,  we  subtracted  0.47  from  the  entries  for  t. 
The  agreement  between  the  observed  values  of  6m,  viz.  6^^  and  the  theoretical 
values,  6fl^,  is  as  good  as  could  be  expected.  This  agreement  appears  to  support 
the  melt-wave  model,  at  least  as  modified  in  this  report,  but  the  theory  should  be 
applied  to  a  number  of  similar  experiments  to  confirm  its  value.  One  remaining 
task  is  that  of  developing  a  theory  for  the  origin  of  tg. 


40 


Acknowledgement 


The  values  reproduced  in  the  Table  above  were  obtained  with  the  assistance 

of  Andrew  Yeoh,  who  help  was  appreciated. 

References 

1.  Woods,  L.C.  “The  current  melt  -wave  model”.  IEEE  Trans,  on  Magnetics, 
Vol  32,  No.  1,  Jan.  1997. 

2.  Woods,  L.C.  lAT  Progress  Report  No.  2,  The  Contact-Spot  Model  of  Transi¬ 
tion,  Dec.  1995. 

3.  Woods,  L.C.  lAT  Progress  Report  No.  8,  Viscous  Effects  on  the  Liquid  Inter¬ 
face  Between  Armature  and  Rail,  Nov.  1996. 

4.  Carslaw,  H.S.  &  Jaeger,  J.C.  Conduction  of  Heat  in  Solids,  (2nd  Ed.),  Oxford 
Uni.  Press,  1959  (see  pp.  268  and  295). 

5.  Rosenthal,  L.  Trans.  Amer.  Soc.  Mech.  Engrs.,  68,  849-66,  1946. 


41 


(Second  CollecUon  of  Progress  Deports  on  Pesearch  into 
Electromagnetic  Launchers 


Investigatiofl  of  Transverse  6triations  on 

Aluminum  Pails 


March  1997 


Professor  L.  C.  Woods  Dr.  Andrew  Y<soh 

University  of  Oxford  InsUhitc  for  Advanced  Technology 

MathemaUcal  InsUtute  University  of  Texas 


Investigation  of  Transverse  Striations  on  Aluminium  Rails 

Reference  Subject  Area  9:  Materials  and  Structures 

Leslie  C.  Woods  and  Andrew  Yeoh 
Institute  for  Advanced  Technology 
4030-2  W.  Braker  Ln.  Suite  200 
Austin,  TX  78759 

Tel:  (512)471-9060;  Fax:  (512)471-9096; 

E)-mail:  leslie.woods@balliol.oxford.ac.uk 


Summary: 

The  ohmically-heated,  liquid  layer  between  the 
rails  and  solid  armature  of  a  railgun  has  a  key  role  in 
the  onset  of  transition,  so  understanding  its  nature 
is  important.  Experiments  in  which  both  the  rails 
and  the  armature  are  made  of  aluminium,  show  pe¬ 
riodic,  transverse  striations  in  the  rail  surface,  with 
a  wave  length  varying  between  0.5  mm  and  about 
5  mm.  Apparently  the  melted  aluminium  at  the 
contact  surface  has  solidified  into  waves  with  their 
crests  at  right  angles  to  the  motion  of  the  armature. 
It  is  conjectured  that  these  waves  are  generated  by 
an  instability  in  the  sheared  fluid  layer  separating 
the  rail  and  armature,  at  a  point  just  to  the  rear  of 
the  armature.  While  layers  of  uniform  vorticity  ad¬ 
jacent  to  a  rigid  wall  are  known  to  be  linearly  stable, 
there  is  a  non-linear,  secondary  instability  capable 
of  generating  the  observed  waves.  This  instability 
can  also  cause  ‘filamentation’,  namely  the  ejection 
of  thin  filaments  of  vortical  fluid  from  regions  of 
high  curvature.  That  the  observed  crests  are  rather 
steep,  appears  to  be  due  to  this  phenomenon. 

1.  Introduction 

During  an  electromagnetic  launch,  bore  deposits 
are  usually  left  behind  by  the  armature.  In  the  case 
of  an  aluminium  armature,  a  primarily  aluminium 
deposit  is  found  on  the  recovered  rail  conductors,  re¬ 
gardless  of  the  rail  conductor  material.  The  general 
nature  of  these  deposits  have  been  previously  re¬ 
ported  [1,2];  typically  one  sees  a  smooth-appearing 
layer  (though  microscopically  rough)  that  averages 
20  /xm  thick.  It  should  be  noted,  however,  that  these 
results  originated  from  either  a  monolithic  or  a  lam¬ 
inated  armature,  with  equal  spacing  between  the 
conducting  and  insulating  laminates.  More  partic¬ 
ularly,  the  armature  conductor(s)  spanned  the  rail 
height  almost  completely  in  the  above  experiments. 

The  observations  and  resulting  theory  presented 
in  this  paper  arc  due  to  not  using  the  entire  rail 
conductor  height  during  the  launch.  Specifically, 


only  40%^  of  the  rail  conductor  height  about  the 
longitudinal  centerline  of  the  rail  was  in  contact  with 
the  armature.  The  initial  reason  for  this  geometry 
was  to  enable  a  12.7  mm  x  25.4  mm  armature  to 
be  launched  from  a  25.4  mm  square  bore  in  order 
to  obtain  higher  exit  velocities.  Higher  velocities 
were  attained,  but  then  transverse  striations  on  the 
rails,  lying  orthogonal  to  the  velocity  vector,  were 
also  observed.  These  striations  occurred  when  an 
Al  7075-T6  armature  was  launched  on  A1  6061-T6 
rails,  and  were  not  found  with  other  rail  conductor 
materials. 

2.  Experimental  Arrangement  and  Results 


Fig.  1.  12.7  mm  high  x  25.4  mm  wide  armature  used  in 
the  first  experiment. 

These  experiments  were  conducted  at  the  Insti¬ 
tute  for  Advanced  Technology  Laboratory  located 
in  Leander,  Texas.  The  testbed  used  was  the  re¬ 
built  OAT  or  Okaloosa  Armature  Tester  (renamed 
LeOAT  for  Leander-Okaloosa  Armature  Tester) 
coupled  with  a  3  MJ  power  supply  from  the  for¬ 
mer  Thunderbolt  program  [3],  In  the  first  experi¬ 
ment,  the  testbed  was  cored  to  produce  a  rectangu¬ 
lar  bore  measuring  12.7  mm  high  x  25.4  mm  wide 

^  Althoiiglj  the  bore  size  is  25.4  mm  x  25.4  mm,  the  rail 
height  is  actually  31.8  mm,  of  which  25.4  mm  is  available  to 
the  armature. 
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(the  larger  being  the  rail-to-rail  dimension).  Due  to 
testbed  constraints,  the  rail  height  was  maintained 
at  31.8  mm,  of  which  12.7  nun  was  in  contact  with 
the  armature.  The  armature  used  for  tliis  ex[)eri- 
rnent  is  depicted  in  Fig.  L.  It  had  six  adjacent  2 
mm  thick  .41  7075-T6  laminates  which  produced  a 
12  mm  stack  held  together  by  an  A1  6061-TG  pin 
and  a  polycarbonate  bore  rider.  The  total  mass  of 
this  package  was  19.00  grams. 

For  the  second  experiment,  a  25.4  mm  square  bore 
w^as  used.  However,  the  six  Al  7075-T6  laminates 
were  arranged  in  a  manner  similar  to  the  first  ex¬ 
periment  and  is  illustrated  in  Fig.  2.  The  initial 
mass  of  this  armature  was  26.70  grams. 

ARMATURE  SIDE  VIEW 


Polycarbonate 


Al  7075T6 


Polycarbonate 


Fig.  2.  Schematic  side  view  of  the  25.4  mm  x  25.4  mm 
armature  used  in  the  second  experiment. 

In  both  cases,  Al  6061-T6  rails  were  used  and 
a  driving  current  with  a  450  kA  peak  propelled 
the  armatures.  Because  of  the  initially  smaller 
mass,  the  first  armature  attained  an  exit  velocity 
of  ~2.0  km/s,  while  the  second  armature  reached 
~1.8  km/s.  Transition  occurred  almost  immedi¬ 
ately  (several  centimetres  beyond  x  —  0  cm)  and 
transverse  striations  were  observed  soon  after  the 
location  of  armature  transition  on  the  recovered 
rails.  Figs.  3  and  4  are  photographs  depicting  the 
transverse  striations  formed  by  the  armature  de¬ 
posit  from  each  of  the  experiments. 


Bore 

rider 


25.4  mm 


Fig.  4.  (a)  Fine  transverse  striations  starting  from  the 
x  =  l5  cm  location  and  (b)  coarse  transverse  striations 
about  the  x=30  ern  location.  The  armature  veloci¬ 
ties  were  ~825  ni/s  and  -^1140  rn/s  respectively.  Pho¬ 
tographs  are  from  the  lu'gative  .41  6061-T.G  rail  of  the 
second  exp(‘nm(nit 


It  should  be  noted  that  the  striations  are  by  no 
means  symmetrical,  i.e.  they  do  not  appear  on  both 
rails  simultaneously.  Closer  inspection  of  the  stri¬ 
ated  deposit  revealed  the  presence  of  spherical  cav¬ 
ities  within  the  layer — an  indication  that  the  alu¬ 
minium  fluid  flow  prior  to  solidification  was  turbu¬ 
lent  in  nature.  Fig.  5  is  a  low-magnification  micro¬ 
graph  depicting  a  dust, or  of  spherical  cavities. 

An  analysis  of  these  striations  show  that  the  peak- 
to-pcak  distance  generally  increases  with  the  veloc¬ 
ity  of  the  armature.  The  wavelengths  range  from 
0.5  mm  to  greater  than  5  mm,  but  not  necessarily 
in  that  order,  i.e.  a  1  nim  wavelength  can  be  found 
in  a  region  displaying  [)iimarily  2.5  mm  striations. 
The  tiii(:kne.ss  of  the  de[)osit.  was  al.so  mciisured  to 
vary  from  sev(‘ral  microns  thick  at  the  valleys  of 
the  striations  to  greater  tiian  200/un  thick  at  the 
peaks*  Du(^  to  the;  ext.rcnu^  roughiuvss  of  the  stri¬ 
ated  surfa(’(\  a  more'  accurat(‘  nu'asurement  was  not 
possible. 


Fig.  3.  Ti  ansverso  striations  about  the  x=20  cm  location;  ar¬ 
mature  velocity  was  850  m/s  at  this  location.  Photograph 
from  tlie  positive  AL  6061-TG  rail  of  the  first  experiment. 


“The  rail  liaii  a  iioitiinal  tiiickness  ol  f».35  iiini,  witli  a 
loleranc.e  of  i  25  /on 
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Fig.  5.  Clusters  of  spherical  cavities  within  the  deposit  layer 
suggests  turbulent  mixing  during  the  moment  of  solidification 


3.  The  Solidification  Time 

It  is  a  reasonable  assumption  that  the  waves  ob¬ 
served  on  the  aluminium  rail  are  generated  by  an 
instability  in  the  highly  sheared  fluid  layer  leaving 
the  interface  at  the  rear  of  the  armature.  Vortex 
sheets  are  well-known  to  be  unstable,  a  phenomenon 
termed  the  ‘Kelvin-Helmholtz’  instability.  Classic 
examples  are  observed  in  the  wing-tip  streamers 
generated  by  high-flying  aircraft  in  supercooled,  at¬ 
mospheric  conditions  and  in  the  generation  of  waves 
at  the  free  surface  of  a  liquid,  over  which  a  gas 
is  flowing.  But  these  flows  are  linearly  unstable, 
whereas  sheared  flow  past  a  rigid  wall  is  stable  to 
linear  perturbations.  There  is  a  more  subtle  desta¬ 
bilising  mechanism,  as  we  shall  describe  shortly. 

Viscosity  very  slowly  eliminates  the  fluid  shear, 
but  as  we  shall  show,  in  the  conditions  of  the  I  AT 
experiment,  long  before  this  can  happen,  heat  is  lost 
to  the  rail  and  the  liquid  aluminium  solidifies.  To 
verify  these  statements,  we  need  the  time  tt,  =  h} /u 
for  the  kinematic  viscosity  i/  to  remove  the  shear 
in  a  layer  of  thickness  h.  With  the  typical  value, 
u  =  5.7x10“*^  m^s”^  for  molten  aluminium  and  the 
measured  layer  thickness  of  h  =  50  /xm,  we  find  that 
^  4.4  ms.  We  shall  find  below  that  the  time  for 
the  solidification  of  the  layer  is  much  smaller  than 
this,  so  we  can  ignore  the  effects  of  viscosity  for  the 


Fig.  G,  Solidification  Distance,  i) 


liquid  layer. 

There  remains  the  possibility  that  the  layer  solid¬ 
ifies  before  the  non-linear  instability  that  we  shall 
discuss  below,  has  had  time  to  develop  a  substantial 
perturbation  on  the  fluid  surface.  So  we  need  an  es¬ 
timate  for  the  solidification  time  is  for  the  molten 
liquid  to  give  up  its  heat  to  the  the  solid  vail  surface. 

In  Fig.  6  we  show  an  armature  extending  over  a 
contact  distance  a,  measured  in  the  direction  of  mo¬ 
tion.  The  mixture  remains  liquid  on  the  rail  surface, 
until  a  distance  (a  -f  D)  measured  from  the  leading 
edge  of  the  armature.  We  term  D  the  solidification 
distance.  Denote  by  Tm^To  the  melting  temperature 
and  the  laboratory  temperature,  and  let  Ar  denote 
the  rail  thickness.  The  surface  t/  =  0  separates  the 
liquid  phase  from  the  solid  phase.  Ignoring  ohmic 
heating  in  the  rail,  we  find  that  its  temperature  is 
approximately 

T-To={l-y/Xr){T,n-To)  (-Ar  <  y  <  0). 

(1) 

Let  the  rail  be  moving  with  a  velocity  V  relative 
to  the  armature  in  a  direction  along  the  OX-axis. 
Also  let  Kr  denote  the  thermal  conductivity,  p  the 
density  and  Lm  the  enthalpy  of  fusion.  In  the  range 
a  <  X  <  00,  the  aluminium  layer  is  not  ohmically 
heated  and  viscous  heating  is  negligible.  Hence,  in 
steady  conditions,  the  energy  equation  becomes 

q2t 

0  =  /Cr-^  4-  pLm{\V)8[x*)  (x*  =  a  -h  D), 

where  (5(x)  is  the  delta  function,  the  liquid/solid  in¬ 
terface  in  the  layer  being  at  x  =  a  +  Z?  with  an 
average  velocity  of  relative  to  the  liquid  layer. 
First  we  integrate  across  the  liquid  layer  to  get 

At  y  =  0,  =  -(r„.  -  To)/Ar;  at  y  =  /i,  «  0, 

since  we  shall  ignore  losses  to  the  air.  Hence 

Y{T^-To)  =  \hpLmV6{x*). 

Next  we  integrate  from  x  =  atox  =  a-l-D.  This 
gives  the  solidification  distance 

where  we  have  used  the  values  (in  MKS  units) 
Kr  =  251,  p  =  2410,  L™  =  3.97x10®,  T„,-Tq  =  635. 
In  the  experiment,  V  =  850  m/s,  Ar  =  6.35  mm, 
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h  =  50 /xm  and  therefore  from  (2),  £)  «  809mm. 
The  solidification  time,  t  =  D/V^  is  therefore  about 
0.95  ms.  This  time  is  to  be  compared  with  the 
growth  time  of  the  non-linear  shearing  instability, 
to  be  described  below. 

4.  Instability  of  Finite  Amplitude  Waves  in 
a  Sheared  Fluid 

The  instability  in  a  sheared  fluid  has  a  long  his¬ 
tory,  that  started  with  Kelvin  and  Helmholtz.  They 
considered  vortex  sheets,  which  work  was  extended 
by  Lord  Rayleigh  [1]  to  sheets  of  finite  thickness. 
Rayleigh’s  unstable  sheets  were  free  to  distort  sinu¬ 
soidally,  maintaining  constant  width.  The  greatest 
growth  rates  are  found  at  wave  lengths.  A,  equal  to 
about  8  times  the  thickness  h  of  the  sheet.  His  lin- 
eai  theory  showed  that  infinite  sheets  adjacent  to 
rigid  walls  were  not  unstable,  at  least  to  small  dis¬ 
placements. 

Consider  the  case  of  a  sheet  of  thickness  /i,  in 
which  the  fluid  velocity  deceases  uniformly  from  a 
value  V  directed  along  the  positive  x-axis,  on  the 
solid  surface  y  =  — /i,  to  zero  on  y  =  0.  The  result¬ 
ing  vorticity  is  u  =  V/h  and  each  fluid  element  spins 
with  an  angular  velocity  ^  The  wail  moves 

with  the  velocity  V,  although  as  we  shall  ignore  vis¬ 
cosity,  this  is  not  an  essential  element.  Above  y  =  0 
the  flow  is  irrotationai.  Rayleigh  adopted  a  theo¬ 
rem  due  to  Helmholtz  that  stated  that  the  effect  of 
an  element  dA  rotating  with  an  aingular  velocity  C, 
is  to  produce,  at  a  point  whose  distance  from  the 
element  is  r,  a  transverse  velocity  g,  such  that 


Fig.  7.  Kelvin-Helmholtz  waves 

It  is  not  difficult  to  show  that  waves  of  length  A 
on  the  surface  y  =  0  will  propagate  with  a  velocity 
c  given  by 

c=i(I-e-.)  (4) 


These  are  called  Kelvin-Helmholtz  waves;  they  have 
neutral  stability.  In  recent  years,  with  the  aid  of 
computers,  this  work  has  been  extended  to  the  non¬ 
linear  regime  [5,6].  For  a  general  survey,  Saffman’s 
monograph  [7]  may  be  consulted. 

The  configuration  of  interest  here  is  that  of  a 
5cmi-infinite  sheet,  extending  from  the  rear  of  the 
^mature  at  x  =  0  to  the  breech  at  (say)  x  ~  oo. 
As  illustrated  in  Fig.  7,  the  K-H  waves  will  prop¬ 
agate  away  from  the  armature,  but  the  presence 
of  the  armature  boundary  changes  the  stability  of 
these  waves.  Consider,  for  example,  the  effect  of  an 
small  increase  in  the  wave  amplitude  at  the  point 
a  of  the  wave  closest  to  the  armature.  A  second- 
order  effect  of  this  modulation  will  be  to  reduce 
the  strength  of  the  vorticity  w,  and  as  this  effect 
is  convected  with  the  local  fluid  velocity,  it  is  ev¬ 
ident  from  (3)  that  the  reduced  spin  of  the  fluid 
elements  immediately  downstream  of  the  point  a 
will  tend  to  increase  the  amplitude  of  the  pertur¬ 
bation  at  a.  In  an  infinite  wave  pattern,  perturbed 
over  all  waves,  this  tendency  would  be  balanced  by 
an  increase  in  the  spin  at  the  trough  of  the  wave 
immediately  upstream  of  a.  Our  conclusion  is  that 
the  presence  of  the  armature — which  suppr^ses  up¬ 
stream  changes — induces  an  instability  in  the  nat¬ 
urally  occurring  K-H  waves  in  the  layer  behind  the 
armature.  We  have  not  attempted  to  quantify  this 
phenomenon,  since  there  are  other  non-linear  effects 
to  take  into  account.  The  flow  pattern  is  rather 
complex.  Just  as  in  the  linear  theory  for  free  vortex 
sheets,  there  should  be  a  value  of  the  ratio  A//i  at 
which  the  growth  rate  of  the  instability  is  a  maxi¬ 
mum  and  so  we  expected  the  observations  to  show 
that  X/h  =  a,  where  a  is  roughly  constant.  And 
certainly  the  longer  waves  occurred  where  the  layer 
was  thickest,  but  no  clear  value  of  a  emerged. 

Pullin  [5]  shows  that  convecting  free  vortical 
boundaxies  are  quickly  distorted  by  differential  mo¬ 
tions.  For  large  enough  wave  amplitude,  the  val¬ 
leys  in  the  wave  pattern  tend  to  move  faster  than 
the  peaks,  a  phenomenon,  one  can  slso  trace  to  the 
influence  of  the  distorted  vorticity  distribution,  as 
indicated  in  Fig.  7.  If  this  process  occurs  in  the  alu¬ 
minium  layer,  we  would  expect  to  find  the  crests  of 
the  waves  to  be  thicker  towards  the  muzzle  end  of 
the  rail.  This  is  observed  to  be  the  case  in  Fig.  4(a), 
although  at  another  station  nearer  the  muzzle,  the 
opposite  wave  inclination  appears. 

Pullen  also  shows  that  the  time  scale  for  the 
distorsion  of  the  primary  wave  is  approximately  the 
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same  as  the  wave’s  period,  viz. 

(jJ 

In  a  typical  railgun  observation,  with  h  = 
50  /xm,  A  =  800  /xm  and  uj  =  V/h  = 

850/(50x10-®)  =  1.7x10^  we  get  tm  =  1.3 /xs, 
which  is  somewhat  less  than  the  soldification  time 
calculated  above,  i.e.  the  waves  have  time  to  de¬ 
velop. 

Finally  we  should  mention  the  phenomenon  of  fil- 
amentation  studied  by  Pullen  et.  ai.  This  is  the  ten¬ 
dency  of  small  bumps  on  surfaces  separating  regions 
of  differing  vorticity  to  steepen  and  form  points  of 
high  curvature,  from  which  spring  thin  filaments  of 
vortical  fluid.  There  is  some  evidence  of  this  with 
the  aluminium  waves,  the  crests  of  which  sometimes 
have  a  rough  ‘Christmas  tree’  edging. 

4.  Conclusions 

It  is  not  easy  to  drawn  firm  quantitative  conclu¬ 
sions  from  the  experiments,  because  of  the  variabil¬ 
ity  of  the  ratio  X/h.  The  thickness  h  changes  over 
an  order  of  magnitude  from  one  part  of  the  rail  to 
another.  That  the  various  non-lineax  instabilities  we 
have  cited  are  involved  in  the  process,  we  have  no 
doubt.  The  evidence  about  the  direction  of  cresting 
is  a  little  confusing,  but  that  filamentation  is  occur¬ 
ring  seems  convincing.  It  seems  unlikely  that  alu¬ 
minium  rails  have  any  practical  value,  except  per¬ 
haps  to  be  a  means  of  studying  the  very  complicated 
liquid  metal  boundary  between  armatures  and  rails. 
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